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ABSTRACT 

We present the first radial velocity measurement of the stellar component of the 
Local Group dwarf galaxy Phoenix, using FORS1 at the VLT UT1 (ANTU) telescope. 
From the spectra of 31 RGB stars, we derive an heliocentric optical radial velocity of 
Phoenix Vq = —52 ± 6 km s" 1 . On the basis of this velocity, and taking into account 
the results of a series of semi-analytical and numerical simulations, we discuss the 
possible association of the HI clouds observed in the Phoenix vicinity. We conclude that 
the characteristics of the HI cloud with heliocentric velocity -23 km s _1 are consistent 
with this gas having been associated with Phoenix in the past, and lost by the galaxy 
after the last event of star formation in the galaxy, about 100 Myr ago. Two possible 
scenarios are discussed: the ejection of the gas by the energy released by the SNe 
produced in that last event of star formation, and a ram-pressure stripping scenario. 
We derive that the kinetic energy necessary to eject the gas is EsNe ~2x 10 51 erg, 
and that the number of SNe necessary to transfer this amount of kinetic energy to the 
gas cloud is ~ 20. This is consistent with the number of SNe expected for the last 
event of star formation in Phoenix, according to the star formation history derived 
by Martmez-Delgado et al. (1999). The drawback of this scenario is the regular 
appearance of the HI cloud and its anisotropic distribution with respect to the stellar 
component. Another possibility is that the HI gas was stripped as a consequence of 
ram-pressure by the intergalactic medium. In our simulations, the structure of the gas 
remains quite smooth as it is stripped from Phoenix, keeping a distribution similar to 
that of the observed HI cloud. Both in the SNe ejection case and in the ram-pressure 
sweeping scenario, the distances and relative velocities imply that the HI cloud is not 
gravitationally bound to Phoenix, since this would require a Phoenix total mass about 
an order of magnitude larger than its total estimated mass. Finally, we discuss the 
possibility that Phoenix may be a bound Milky Way satellite. The minimum required 
mass of the Milky Way for Phoenix to be bound is M M w(< 450kpc) > 1.2 x 10 12 M 
which comfortably fits within most current estimates. 

Subject headings: galaxies: individual (Phoenix), galaxies: dwarf, galaxies:evolution, 
ISM: HI, Local Group 

1. Introduction 

At the low end of the mass and luminosity distributions, dwarf galaxies are the most common 
type of galaxies in the Universe. They are usually divided in two groups: dwarf irregular (dlrr) 
and dwarf spheroidal (dSph) galaxies. In the simplest scheme, the key differences between the 
two groups is that dlrr have large amounts of gas and substantial current star formation activity, 
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while dSph lack gas and recent star formation. However, this picture has been shown to be too 
simplistic, and the distinction between the two classes is not always straightforward. There are 
dlrr galaxies, like LGS 3 (Aparicio, Gallart & Bertelli 1997b; Young &; Lo 1998) and Pegasus 
(Aparicio, Gallart & Bertelli 1997a; Lo, Sargent & Young 1993), with substantial amounts 
of gas but apparently low star formation rate. Other galaxies, like Fornax and Leo I, which 
are traditionally considered dSph, have relatively recent (few hundred Myr old) star formation 
(Stetson, Hesser & Smecker-Hane 1998; Gallart et al. 1999a,b). Finally, there are some interesting 
cases where gas, not centered in the optical component of the galaxy, is or may be associated with 
dSph systems: Sculptor (Carignan et al. 1998), Tucana (Oosterlo, Da Costa & Staveley-Smith 
1996), And III, And V, Leo I and Sextans (Blitz & Robishaw 2000). In the case of Sculptor, the 
kinematic association of the galaxy with two HI clouds located symmetrically about the galaxy 
center has been confirmed by Carignan et al. (1998). 

From the theoretical point of view, the relationship between dlrr and dSph galaxies, and their 
possible evolutionary links have been profusely discussed in the literature (e.g. Einasto et al. 1974; 
Lin & Faber 1983; Dekel & Silk 1986; Silk, Wyse k Shields 1987; Davies & Phillips 1988; Mac 
Low &: Ferrara 1999; Ferrara & Tolstoy 2000). The fact that the same basic relations and overall 
structure are shared by dlrr and dSph (Bingelli 1993) seem to point to a common origin. On the 
other hand, the morphological segregation of dwarfs, with dSph in general close to large galaxies 
and preferentially isolated dlrr (van den Bergh 1994) may suggest some kind of environmentally 
driven evolution, either related to ram pressure sweeping (Lin & Faber 1983), or enhanced star 
formation that rapidly consumed the gas in high density environments (Davies &: Phillips 1988; 
Lacey & Silk 1991). 

Phoenix has the characteristics of a system in which a transition from a dlrr galaxy to a 
dSph galaxy may be currently going on. It shows recent star formation in its central part (500 
to 100 Myr ago), but no HI has been observed centered in the optical galaxy. Nevertheless, 
several clouds at different radial velocities have been detected close to it. Carignan, Demers & 
Cote (1991) detected HI emission at Vq = 56 km s _1 which was clearly separated from a much 
larger scale component at ~ 120 km s _1 associated with the Magellanic Stream. Oosterloo et al. 
(1996) detected another smaller cloud at Vq= -23 km s~ x , situated at 6' SW from the center 
of the galaxy. VLA observations by Young & Lo (1997) confirmed the existence of this last 
cloud close to the optical galaxy, forming a curved structure that wraps around its SW border. 
Subsequently, St-Germain et al (1999) reobserved a 2deg x2deg area around Phoenix with the 
Australia Telescope Compact Array, remapping all the previously observed components, plus one 
at 7 km s _1 which is likely of Galactic origin. After a thorough discussion of all available evidence, 
these autors conclude that, due to its compact shape and internal velocity gradient, the amount 
of HI mass (~ 10 5 M ) derived for the cloud if at the distance of the galaxy, and its location, 
almost overlapping the SW outskirts of Phoenix, the -23 km s _1 component is the one most likely 
associated with Phoenix. It is interesting that the spatial variation of the young population with 
age in Phoenix is consistent with a possible self-propagation of the star formation over its disk 
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from East to West (Martinez-Delgado et al. 1999; see also Held, Savianne & Momany 1999). 
Therefore, if the -23 km s _1 component would prove to be associated with Phoenix, it would 
provide tantalizing evidence that a burst of star formation may blow-out the gas from a dwarf 
galaxy. However, without a measurement of the radial velocity of the Phoenix stellar component, 
it is not possible to stablish a reliable connection between the stars and the HI gas seen near 
Phoenix in projection. 

We present here the first measurement of the radial velocity of the stellar component in 
Phoenix, and discuss the possible association between the galaxy and the gas clouds in its 
neighbourhood in the sky. This paper is organized as follows: in Section [2| we present the 
spectroscopic observations performed with FORS1 at VLT-ANTU. The derivation of the radial 
velocities of the Phoenix red giant branch (RGB) stars is described in ||. Finally, the implications 
of the derived Phoenix radial velocity regarding the possible association to the dlrr of the different 
HI clouds detected toward it is discussed in Section ||, based on the results of semi-analytical and 
numerical simulations of SNe explosions and ram-pressure stripping scenarios. Our conclusions 
are summarized in Section [|. 

2. Observations 

Spectroscopy of a sample of candidate Phoenix stars and several velocity standard stars was 
obtained with FORS1 at VLT-ANTU on the night of September 14-15, 1999. The observing 
conditions were apparently photometric, but the seeing was mediocre, about ~ 1-1.5" the whole 
night with some improvement only at the end of the night, down to 0.7". 

We knew from past experience that, even for red giant stars, the best spectral range to 
observe them for radial velocity measurements purposes is between about 450-530 nm (Vogt et 
al. 1995). In this region, even for relatively metal-poor stars, we obtain the highest density of 
spectral lines, adequate continuum stellar flux, and low sky emission. We extended the spectral 
coverage some more to the red to include some sky lines, in particular the [01] 5577 line, to set the 
zero-point of the wavelenght callibration. In FORS1, the grism that covers that spectral region 
is GRIS600B. Its spectral resolution is R=815, which, with a resolution of 50A/mm and a CCD 
pixel size of 24^m, gave a dispersion of 1.2A/pixel. 

We observed in multislit mode (MOS) to obtain spectroscopy of a sample of stars in the field 
of Phoenix (Figure |l|). The usable field of view of 6'. 8 x 4' basically covered the main body of the 
galaxy. In two MOS setups (see Table |l| and Figures ^ and ||), we were able to employ 14 and 
16 slitlets on stars that had a very high probability of membership as red giants in the galaxy, 
as selected from the Martmez-Delgado et al (1999) photometry. A second star happened to be 
aligned with the selected target star in two of the slitlets of the second set. Figure |I| shows a 
Phoenix CMD with the observed stars marked. 

We prepared approximate MOS slitlet configurations using the FIMS VLT software and 
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the LCO image of Phoenix used in Martinez-Delgado et al. (1999). This image did not allow 
positioning the slitlets with enough precision to actually perform the observations, but it did 
allow to select optimal configurations of the 19 slitlet set to secure maximum occupancy by the 
brightest red giants. The night before our observations, Dr. Iovino (who was observing before us 
with FORS1) kindly took a 120 sec V FORS1 image of Phoenix for us, with which we were able 
to refine the slitlet positioning in such a way that no further adjustements were necessary during 
the target adquisition at the telescope. 

Four velocity standards of spectral type K0-K3 were observed a total of 7 times through the 
night. Because the exposure times necessary to obtain good S/N spectra of the standard stars were 
short to provide a good measurement of the position of the 5577 01 line, we took, immediately 
after each standard star measurement, a companion spectra of the sky, with exposure times of 
300-600 sec. These companion spectra proved indeed very useful to obtain reliable corrections of 
the zero point of our wavelenght calibration. Table [A] identifies the standard stars observed, and 
lists the velocity measurements. In addition, we took a 300 sec exposure of the twilight sky to be 
used as an extra, high S/N velocity template. 

Arcs and screen flatfields were taken in the morning after our observations in each MOS setup. 
For practical considerations, the obtention of callibration arcs during the night is not supported 
with FORS1. This is because the whole operation of recording a wavelength calibration spectrum 
at night requires approximately 30 minutes. Because it involves a reconfiguration of the telescope 
and of the instrument, it does not offer any significant advantage, compared to next morning 
daytime calibrations, in terms of stability of the configuration with respect to the one used for the 
observation of the scientific target. Therefore, the use of arc lamp spectra is restricted to defining 
the dispersion relation while zero point alignment is more accurately achieved using the night sky 
lines as references. 

3. Velocity measurements 

3.1. Basic reduction procedures 

Because this is our first paper describing radial velocity measurements using FORS1 at the 
VLT, we will describe our reduction techniques in some detail. This section may be skipped by 
any readers uninterested in the precise reduction procedure. We mainly used routines from the 
EXPORT version V2.ll of Sun/IRAF, with occasional use of MIDAS routines. 

1) The data were taken in single amplifier mode and high gain (0.68 e~/ADU and 5.03 e~ 
read-out noise). We used IRAF to remove the bias level from the images, by using a set of bias 
frames taken in the same mode, plus the overscan region of each image. 

2) Screen flatfields were obtained in each slitlet configuration using the two available sets (blue 
and red filtered) of halogen flatfield lamps. Four flatfields were averaged for each configuration 
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using IRAF, and subsequently normalized using NORM/MOS whithin the MIDAS package. This 
routine averages the rows separately for each slitlet, smooths with a median filter (of 10 pixels 
here) and divides each row of the slitlet by the filtered average. The division of each object frame 
by the corresponding flatfield was performed using IRAF. 

3) Three 1500 sec exposures were obtained in one of the slitlets configurations while two sets 
(1500 + 900 sec and 2x1800 sec) were obtained in the second. No noticeable shifts were detected 
in either the spatial or the dispersion direction, and the spectra in each set were averaged to 
increase S/N and eliminate cosmic rays. 

4) Individual 2-D spectra were cut from the FORS1 image. From these, one-dimensional 
spectra were extracted using the noao.twodspec.apextract.apall package in IRAF. At this point we 
also substracted the sky background by fitting it across the dispersion at either sides of the star 
spectrum. The substracted sky spectrum was kept in one of the beams of the multispec format in 
order to be eventually used at a later phase to refine the wavelenght callibration using sky lines. 

5) He + HgCd lamp callibration spectra in each slitlet configuration were obtained the morning 
after the Phoenix and standard stars' observations, and were extracted, except for the background 
substraction, with identical parameters as the object spectra. We used noao. onedspec. identify to 
identify the lines in each callibration spectrum. Order 4 Chebyschev polynomia were typically 
fitted to the position of 11-14 spectral lines and the derived transformations were applied to the 
program objects using the noao. onedspec. dispcor package. 

6) The wavelenght callibration zero-point was checked and refined using the [OI] line at 
5577.34 A. This line is present with high S/N in all Phoenix long exposure spectra, but it is too 
faint in the radial velocity standard star spectra to provide a reliable measure of its position. For 
them, we used the longer exposure spectra obtained after shifting slightly the pointing of the 
telescope to remove the standard star from the slitlet (see Table |]). The center of the 5577 line 
was obtained by fitting a gaussian profile using noao. onedspec. splot. The fact that there is some 
correlation between the position of the telescope and the amount and sign of the shift (see Table |3|) 
gives us confidence that we are indeed measuring an actual effect of telescope or instrument 
flexures on the spectra. The calculated 8X shifts, typically within ± 0.4 A(Table |3|) were applied 
to both Phoenix and standard stars' images, by modifying the CRVAL1 header parameter, which 
provides the starting wavelength of each spectra. 

7) The spectra were continuum-substracted using noao. onedspec. continuum. Finally, for the 
Phoenix spectra, we removed all pixels with values > 40 or < -40 by replacing them by 0.0 using 
images. imutil.imreplace to get rid of any sky residuals and remaining cosmic rays. 
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3.2. The Master Radial Velocity Template 

We combined the standard star's observations plus the twilight solar spectra into a single, 
high S/N, master template spectrum. To do this, we first transformed these spectra to a log A scale 
using noao.onedspec.dispcor. We then used the cross-correlation task noao.rv.fxcor to determine 
the relative velocities between each standard star spectrum and an arbitrarily chosen reference 
spectrum (HD 1 76047^3). The relative velocities V re i were transformed to <5(logA) increments as 

<5(logA) = -V rc i x loge/c 

and added to the CRVAL1 header parameter, which, in a logarithmic wavelenght scale, provides 
the logarithm of the starting spectrum wavelenght. Then, we combined the individual spectra 
into a single, high S/N template using noao.onedspec.scombine. By construction, the "heliocentric 
velocity" of this master template is that of HD 176047, i.e. -42.8 km s _1 . 

We then used the master template to obtain the observed radial velocity of the velocity 
standards. From the differences between the observed and true heliocentric velocity, we can i) 
calculate an offset between the radial velocity template and the overall velocity system defined 
by the standards, which will be applied to the Phoenix velocity measurements to transform them 
to a true heliocentric system, and ii) estimate the mean error with which we are determining 
the velocities of the standard stars. Because of the higher S/N of the standard star spectra, this 
will be a lower limit of the error in the Phoenix velocity determinations. Only the standards for 
which we had a companion long exposure sky spectra to correct the zero-point of the wavelenght 
callibration have been used for this purpose. Table [A] lists the results of this comparison. The a v 
values have been calculated as described in Section |3.3| . The cr-weighted mean difference between 
the true and observed heliocentric velocity is Vheiio — Vobs = 1-2 with dispersion oy=1.7 km s _1 . 

3.3. Velocity of Phoenix giants 

The radial velocities of the Phoenix giants were determined by cross-correlating each star 
with the master template described in section |3.2| using the IRAF task noao.rv.fxcor. Figure [5] 
shows the template spectra, an example of a typical Phoenix spectrum, and the corresponding 
cross-correlation function. 

The a v errors were determined, in the same way as described in Vogt et al. (1995) using six 
repeated measurements of two standard plus two repeated measurements of 6 Phoenix stars, listed 
in Table |I[ We assumed that the velocity errors can be written as a v = a/(l + R), where a is a 
constant (Tonry & Davis 1979). We then defined the x 2 statistic for the repeat measurements as 
XX 1 + Ri) 2 ( v i— < v >) 2 /a 2 where < v > refers to the mean velocity of the star corresponding 
to observation i. We then calculate the a value that gives a value of x 2 corresponding to a 50% 
probability that we would exceed that value of x 2 by chance (we will denote this value as xio)- 
With this definition, the adopted errors are equivalent to la Gaussian errors. For this set of 18 
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observations of 8 stars with 10 degrees of freedom, we have xio = 0.934, from which ct=455. 

Table |5| lists the measured individual heliocentric velocities for the observed Phoenix stars 
with the corresponding R and a v values. Only stars with well defined cross-correlation peaks have 
listed velocities. Some stars in Phoenix set#2 have repeated measurements. Two of the observed 
stars turned out not to be RGB stars. Their spectra is briefly discussed in Appendix A. 

Figure || shows that there is no noticeable trend in the measured velocities as a function of 
V magnitude, (V — I) color index, position angle PA, or galactocentric distance r in Phoenix. 
The ovweighted average of these velocities is -57 db 6 km s _1 , the same as the median velocity. 
The mean value of the Phoenix velocity remains stable within the errors if we eliminate from the 
average calculation stars with velocities 3a above and below the mean. With this selection, we 
obtain Vphoenix-obs = — 53 dz 6 km s _1 . After correcting for the difference between the true and 
measured heliocentric velocity of the standard stars, we conclude that the true radial velocity of 
the stellar component in Phoenix is Vq = —52 dz 6 km s _1 , which we will adopt in the remaining 
of this paper. 



4. Discussion. 

4.1. The stellar radial velocity of Phoenix and its association with the HI clouds. 

In their paper about the HI in the field of Phoenix, St-Germain at al (1999; see their Figure 4) 
discuss the identification of four components at heliocentric velocities -23, 7, 59 and 140 km s _1 . 
The components at 7 and 140 km s _1 are most likely of Galactic and Magellanic Stream origin 
respectively. The nature of the other two components is unclear but St-Germain et al. (1999) 
conclude that the component at -23 km s -1 is most likely associated with Phoenix because of i) 
its compact structure and partial overlap with the Phoenix stellar component, ii) the amount of 
HI mass inferred for the cloud at the Phoenix distance, which is similar to the mass of other dSph 
or dSph/dlrr galaxies, hi) its velocity structure, unlike that of HVC, which could indicate either 
rotation or ejection from Phoenix and iv) the location of the cloud west of Phoenix, consistent 
with the propagation from East to West of the most recent star formation event in the galaxy, as 
discussed in Martinez-Delgado et al. (1999). Our determination of the stellar radial velocity of 
Phoenix Vq = —52 ± 6 km s _1 tends to support this particular association, if any, between HI in 
the field of Phoenix and the dwarf galaxy itself. The relatively large velocity difference of gas and 
stars, however, means that the association is not unquestionable and merits some discussion, as 
does the particular morphology and location of the HI cloud. 

In this section we will present the results of a series of semi-analytical and numerical 
simulations, under the assumption that the HI cloud at -23 km s _1 was once part of Phoenix and 
participated in the last star formation event observed in the galaxy about 100 Myr ago, which 
is the age of the younger HeB stars, found in the west region of Phoenix. We will consider two 
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possible mechanisms that may have stripped Phoenix of its gas: gas ejection by SNe explosions 
and ram-pressure stripping. These models allow us to set limits to the necessary conditions for this 
stripping to have taken place under a given mechanism, and we find out that both mechanisms 
may indeed have been responsible of the observed situation. Similar calculations under the 
assumption that the HI cloud at +59 km s" 1 was once associated to Phoenix lead to unlikely 
boundary conditions (very large number of SNe required or large density of the intergalactic 
medium, IGM), which makes us conclude that this one is a chance superposition. 

4-1.1. Gas ejected by SNe explosions 

The possibility that an event of star formation in a small galaxy can actually blow out all or 
part of its gas, has been profusely discussed in the literature from the theoretical point of view 
(e.g. Mathews & Baker 1971; Larson 1974; Dekel & Silk 1986; Mac Low & Ferrara 1999). This gas 
removal could eventually halt subsequent star formation and explain the low metallicity observed 
in dwarf galaxies and the relationship between their total mass and metallicity. Also, the HI 
mappings of some dwarf galaxies (Puche & Westpfahl 1993) have shown that, in the smallest dwarf 
irregular galaxies, one large, slowly expanding shell usually dominates the interstellar medium. 

In the case of Phoenix, we estimated the energy that should have been injected to the 
interstellar medium by SNe explosions if these were responsible for the ejection of the gas. We 
assumed that the HI cloud was ejected 100 Myr ago and projected in the opposite direction with 
respect to the observer (as indicated by the respective observed velocities -52 km s _1 and -23 
km s _1 of Phoenix and HI cloud) and along the trajectory of Phoenix, in such a way that both 
observed velocities are affected by the same projection angle. We solved the equations of motion, 
considering the only effect of Phoenix's gravity force and imposing the projected distance between 
Phoenix and the HI cloud (d ~ 0.65 kpc) and their observed radial velocities (-52 km s _1 and 
-23 km s _1 , respectively) as boundary conditions. This calculation allows us to simultaneously 
estimate the energy necessary to eject the gas, E$Ne ~ 2 x 10 51 erg and the most likely angle 
between the line of sight and the line connecting Phoenix and the HI cloud, a ~ 10 degrees. 

To estimate the number of SNe necessary to inject this energy into the interstellar medium, 
we should estimate what fraction of the energy of a SNe is actually transferred as kinetic energy 
to the gas. Larson (1974) showed that the thermal energy available per SNe for driving a galactic 
wind in cases where SNe remnant cooling is important (as would be the case for dwarf galaxies) 
is of the order of 10 50 erg, i.e. about 10% the total energy output of the SNe. In the case of 
Phoenix, therefore, about 20 SNe would have been necessary to produce the inferred energy. By 
comparing the number of observed 100 Myr old stars in the main sequence of Phoenix with the 
number predicted by synthetic CMDs, Martmez-Delgado et al. (1999) calculate that the total 
mass involved in the burst of star formation that produced the younger stellar association in the 
western part of the galaxy is 1.8 x 10 4 M Q . Using the initial mass function of Kroupa, Tout Sz 
Gilmore (1993), we derive that ~ 50 stars more massive that 1OM should have been produced 
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during that event, providing a potential number of SNe large enough to produce the necessary 
energy to have blown away the gas. 

The most likely estimated angle a ~ 10 degrees between the line of sight and the line 
connecting Phoenix and the HI cloud means that both objects are in almost radial trajectories 
with respect to the line of sight. Therefore, their total velocities are very similar to the projected 
component, and the actual distance between Phoenix and the HI cloud is about 3.5 Kpc. 
Considering this distance and the relative velocities of the HI cloud and Phoenix, a Phoenix 
total mass Mph oe ^ tot ~ 4 x 1O 8 M would be necessary if the HI cloud was to be gravitationally 
bound to Phoenix. This is about one order of magnitude larger than the total mass estimated for 
Phoenix, ~ 3.3 x 1O 7 M (Mateo 1998), which would imply that (M/L) Phoe ~ 4OOM /L Q , and 
indicates that indeed the amount of energy produced by the last burst of star formation would 
have been large enough to blow away the small amount of gas remaining in the galaxy after that 
particular episode of star formation. 

This model of SNe driven winds, even though very convincing from the energetic point if 
view, has the drawback of the regular appearance of the HI cloud, and the highly anisotropic 
geometry of the gas-stars configuration, with all the gas out in one direction in a single isolated 
cloud. These are not a problem for the ram-pressure hypothesis discussed in the next section. 

4-1.2. Ram-pressure stripping 

The gas of the IGM produces a ram-pressure on the gas of an object (satellite) moving 
through it 

Pram OC PIGMV (1) 

where pigm is the gas density of the IGM and v is the velocity of the satellite (Gunn & Gott 
1972) relative to the IGM. When the pressure is strong enough, the gas of the satellite is slowed 
down and stripped from the satellite. We have simulated a gas rich dwarf galaxy moving through 
a gaseous environment, using a semi-analytical approach to model the ram-pressure stripping (see 
Gomez-Flechoso 2000 for details of the numerical models). We have analyzed the conditions for 
the HI cloud to be stripped by the IGM. We assumed an ambient medium of uniform density 
and at rest with respect to the Galaxy, a peak of the gas density of the HI cloud of 4 x 10 19 
atoms/cm 2 and the total luminosity of Phoenix of Lb = 9.5 x 10 5 L Q (St. Germain et al. 1999). 
Various M/L ratios have been considered for the dwarf galaxy in our simulations. As in the 
previous section, we reproduce the observational velocities and positions of Phoenix and the HI 
cloud, assuming that (i) the gas has been stripped along the trajectory, and (ii) the ram-pressure 
started 100 Myr ago. In this case, the forces involved in the trajectory of the cloud are the 
gravitational force of the dwarf galaxy and the ram-pressure stripping force due to the IGM. The 
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latter is produced by the ram-pressure given by eq. (|l|), taking into account that the velocity 
v in the equation is the velocity relative to the IGM, which is the Galactocentric velocity. The 
Galactocentric radial components of the velocity of Phoenix and the HI cloud are V^p hoe = —144 
km s _1 , and V^^j = —115 km s _1 , respectively. As in the calculation relative to the SNe in the 
previous section, in order to reproduce the projected distance between Phoenix and the HI cloud 
and their Galactocentric radial velocities, we estimate the angle, a, between the line of sight and 
the trajectory of the satellite, and the gas density of the environment, pigm (assuming it is nearly 
constant). The best model corresponds to a ~ 21 degrees, roughly independent of the M/L ratio 
of Phoenix, and p IGM ~ 2.0 - 3.0 x 10~ 4 atoms/cm 3 , for (M/L) Phoenix = 2-20 M /L . This 
value of a implies that: (i) the orbit of Phoenix is quite eccentric, (ii) the true distance between 
Phoenix and the HI cloud is about 1.8 kpc, and (iii) the total velocity of Phoenix and the HI 
cloud are Vp^ oe — —154 km s _1 and V§f — —123 km s -1 . The value of the gas density of the 
environment, piGM = 2.0 — 3.0 x 10 -4 atoms/cm 3 , obtained in this case is consistent with some 
estimates of the hot gas density of a Local Group halo (Suto et al. 1996; but see also Murali 2000; 
Mulchaey 2000 and references therein). 

In Figure we show the results of the simulation: the gray levels refer to the predicted 
velocity of the gas (darker corresponds to more negative values of the velocity) and the white 
contours correspond to heliocentric radial velocities -27, -24, -21 and -18 km s _1 . The black 
contours correspond to the predicted gas density (3.5, 3.0, 2.5, 2.0 and 1.5 xl0 19 cm -2 ). The 
structure of the gas in the simulations remains quite smooth as it is stripped from the satellite, 
keeping a distribution similar to the observed HI cloud. Note also that a) a velocity gradient in 
the sense that the most extreme differences in velocity are seen farthest from the galaxy, and b) 
the velocity of the stripped gas is closer to the systemic velocity of the Milky Way than that of 
Phoenix are predicted by our model (see also Blitz & Robishaw 2000), in agreement with the 
observed Phoenix parameters. 

These arguments support the hypothesis of the ram-pressure as responsible of the gas 
stripping in Phoenix. However, it may demand Phoenix to be in its first approach to the Galactic 
Center, which may be possible if Phoenix is not bound to the MW or it is in a very long period 
orbit. Given its almost radial orbit, near the pericenter Phoenix must cross zones with a much 
larger gas density than in its current position, and at much larger velocity with respect to the 
IGM. This would imply a much larger ram-pressure stripping, to which the gaseous content of the 
galaxy could have only survived if it originally was very massive and concentrated, unlike at the 
present time. 

There is at least a couple of theoretical works that consider tidal stripping of a gaseous 
satellite approaching to the Milky Way. They suggest that such a survival in the close approach 
may not be completely unrealistic. Mayer et al. (1999, 2001) show that low surface brightness 
satellites of a Milky Way-like galaxy, in eccentric orbits (as typical in hierarchical models of 
galaxy formation), may lose up to 90% of their stellar mass due to tidal stripping. But, could any 
amount of gas remain still bound to the galaxy after one or more pericentric passages? Lin & 
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Murray (1999) argue that "around the dSph galaxies in the Galactic halo, their ejected unbound 
gas follows their orbits analogous to the Magellanic Stream. Gas density declines/enhances near 
peri/apogalacticon as the stream is stretched/compressed (by tidal forces). At apogalacticon, 
where the density is at a maximum and exposure to the Galactic UV flux is at a minimum, the 
diffuse gas can recombine to form atomic and molecular hydrogen and eventually stars" . To assess 
the plausibility of these different alternatives when ram-pressure is included in the case of Phoenix, 
it would be necessary to perform much more detailed simulations, which exceed the scope of the 
present paper. In any case, gas loss will be maximum at the pericenter, where both tidal stripping 
and ram-pressure forces are larger. 



4.2. Comparison with observations of other dwarf galaxies 

In two cases, DDO 210 and LGS3 for which HI superimposed to the stellar component had 
been observed, Blitz & Robishaw (2000) detect a stronger emission component offset from the 
stellar component. They also find HI possibly associated with Leo I, whose position is also offset 
from the stellar component. Leo I and LGS3, for which an optical velocity is also known, present 
interesting similarities with the Phoenix case. In LGS3, HI cloud and galaxy differ in velocity 
by about 50 km s" 1 , which implies that they are not gravitationally bound. In the case of Leo 
I, the velocity centroid of the emission that seems more closely associated with the galaxy (both 
spatially and kinematically) differs by about 30 km s _1 from the velocity of the stellar component. 
Apart from the magnitude of the difference in velocity between gas and stars, Leo I, LGS3 and 
Phoenix are similar also in a) the morphology of the gas relative to the stars and in b) the time 
of the last events of star formation, which is of the order of a few hundred Myr in all of them 
(Phoenix: Martfnez-Delgado et al. 1999; Leo I: Gallart et al. 1999a, b; LGS3: Aparicio et al. 
1997b). For Leo I and LGS3, Blitz & Robishaw (2000) discuss the possibility that the HI has been 
stripped by the ram-pressure of the hot (Galactic or M31) halo gas. In spite of a distance to the 
Milky Way of almost twice that of Leo I, this seems a plausible explanation also for Phoenix, as 



we have discussed in Section 4.1.2 



4.3. Phoenix as a Milky Way satellite, and the mass of the Galaxy 

The minimum required mass of the Milky Way for Phoenix to be bound (based on the values 
of the Galactocentric velocity of Phoenix and its distance to the Galactic center obtained above) 
is Mmw{< 450kpc) > 1.2 x 10 12 Mq. This value confortably fits within most current estimates of 
the Milky Way mass (see Zaritsky 1999 for a review). For example, Zaritsky et al. 1989 conclude 
that Mmw = (1-3 db 0.2) x 10 12 from timing arguments involving Leo I and the Milky Way; on 
another hand, Anguita, Loyola & Pedreros (2000) calculate, from a recent LMC proper motion 
measurement that Mmw = (3 — 4) x 10 12 out to 50 kpc, if the LMC is to be bound to the 
Milky Way (see also Mendez et al. 1999). We conclude, therefore, that Phoenix is very likely the 
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outermost (known so far) Milky Way satellite. 

5. Summary and conclusions 

We present the first measurement of the radial velocity of the Phoenix stellar component, 
using FORS1 in MOS mode at the VLT UT1 (ANTU) telescope. We have obtained spectra of 31 
RGB stars. Two additional stars for which we obtained spectra turned out to be a Carbon star 
and a bright supergiant star (see Appendix A). From the spectra of the RGB stars we derive an 
heliocentric optical radial velocity of Phoenix V & = —52 ± 6 km s _1 . On the basis of this velocity, 
and taking into account the results of a series of semi-analytical and numerical simulations, we 
discuss the possible association of two of the four HI clouds observed in the Phoenix vicinity, 
those at velocities -23 and +59 km s _1 which are not associated to either the Milky Way or the 
Magellanic Stream (St-Germain et al. 1999; Young & Lo 1997; Oosterloo et al. 1996). 

We conclude that the characteristics (mass, distance, velocity relative to Phoenix) of the HI 
cloud with heliocentric velocity -23 km s _1 are consistent with this gas having been associated 
with Phoenix in the past, until the last event of star formation in the galaxy, about 100 Myr ago. 
Two possible scenarios for the loss of this gas are discussed: its ejection by the energy released 
by the SNe produced in that last event of star formation, and a ram-pressure stripping scenario. 
The assumption that the HI cloud at +59 km s -1 was once associated to Phoenix lead to unlikely 
boundary conditions (very large number of SNe required or large density of the IGM), which 
makes us conclude that this one is a chance superposition. 

We derive that the kinetic energy necessary to eject the gas is E$Ne ~ 2 x 10 51 erg. The 
number of SNe necessary to transfer the inferred amount of kinetic energy to the gas cloud is 
consistent with the total mass involved in the last event of star formation ~ 100 Myr ago by 
Martmez-Delgado et al. (1999). Considering the derived distance between Phoenix and the 
HI cloud and their relative velocities, the Phoenix total mass necessary for the HI cloud to be 
gravitationally bound is about an order of magnitude larger than the total mass estimated for 
Phoenix. This indicates that the amount of energy produced by the last burst of star formation 
would have been large enough to blow away the small amount of gas remaining in the galaxy after 
that particular episode of star formation. The main drawbacks of this SNe driven winds scenario 
is the regular appearance of the HI cloud and its anisotropic distribution with respect to the 
stellar component. Another possibility is that the HI gas observed at -23 km s _1 was stripped as 
a consequence of ram-pressure by the IGM. In our simulations, the structure of the gas remains 
quite smooth as it is stripped from Phoenix, keeping a distribution similar to that of the observed 
HI cloud. 

Since both the SNe and ram-pressure stripping scenarios for gas ejection in Phoenix are 
consistent with the observations, a combined scenario may be also plausible. One possibility is 
that supernovae resulting from a burst of star formation diluted the gas in Phoenix and increased 
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its energy enough for it to be swept away by an already existing intergalactic wind. Another 
possibility is that Phoenix may be crossing regions of variable gas density, and that 500-100 Myr 
ago it ran into a relatively dense clump or stream of intergalactic gas that was responsible for the 
sweeping event that presumably terminated -or temporarily halted- the star formation. 

Finally, we discuss the possibility that Phoenix may be a bound Milky Way satellite. The 
minimum required mass of the Milky Way for Phoenix to be bound is Mmw (< 450kpc) > 1.2 x 10 12 
Mq. This value confortably fits within most current estimates of the Milky Way mass. 

It is a pleasure to thank I. Perez, R. Guzman, E. Hardy and S. Zepf for very useful discussions 
on spectroscopic reductions, to R. Larson and L. Mayer for enjoyable conversations about possible 
scenarios for gas stripping, and for a critical reading of the manuscript. J. Mulchaey and S.G. 
Navarro pointed us to useful references on gas in galaxy groups and spectral classification, 
respectively. We also thank A. Iovino, for taking an image of Phoenix the night before our 
observations which allowed us to produce an accurate MOS setup. Support for this work was 
provided by Chilean CONICYT grant 1990638, NASA grant GO-5350-03-93A and by the IAC 
through grant PB3-94. 

A. A young star and a carbon star spectroscopically confirmed (serendipitously) in 

Phoenix 

One of the stars originally selected as a candidate RGB star, Phoenix set#2-6 and Phoenix 
set#2-10s3 (Table |5|) turned out to be other types of stars. 

Phoenix set#2-6 (Figure |||a) has the spectra typical of a Carbon star (see for example 
Cannon, Niss & Norgaard-Nielsen 1980 for a sample spectra and discussion of the features present 
in a C star spectrum). It is one of the two C stars previously confirmed spectroscopically by Da 
Costa (1993). It is the brightest red star in our sample, located above the TRGB (identified in the 
CMD in Figure || with a solid square). 

Phoenix set#2-10s3 (Figure |s|h>) was not in our original candidate list. It was included in 
the extreme of one of the slitlets, and turned out to be a young star. Its spectra clearly shows 
the series of Balmer lines in absorption, and at this low resolution, low S/N, we can classify it 
in the range of spectral types B5-A5 (Jacoby, Hunter & Christian 1984). Using (m-M)o=23.0 
and Ay =0.06 for Phoenix (Martmez-Delgado et al. 1999), the absolute magnitude for this star, if 
it belongs to Phoenix, is My=-2.83 (Table |J), which would correspond to a bright giant star of 
luminosity class II, and probably early A spectral type (Landolt-Bornstein 1982). The position of 
this star in the CMD is consistent with a blue-loop star with Z=0.001 and age 100 Myr, formed in 
the last event of star formation in Phoenix. 
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Fig. 1. — Identification chart for the stars measured in Phoenix. The image is a 120 sec. V exposure 
at the VLT. North is down and East it to the left. 

Fig. 2. — MOS setup configuration for Phoenix set#l (Table [l]) . FORS1 reference stars are circled. 
West is up and North is to the left. 

Fig. 3. — MOS setup configuration for Phoenix set#2 (Table [l]). FORS1 reference stars are also 
marked. North is up and East is to the left 
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Fig. 4. — CMD of Phoenix with the measured stars near the tip of the red giant branch identified 
as crosses. The C star and the bright supergiant (late B- early A spectral type) star (see Appendix 
A) are identified with a filled square and a filled triangle respectively. 
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Fig. 5. — Top: Extracted spectrum of a typical Phoenix giant (Phoenix set#l-14 in Table |5|). 
Middle: Spectrum of the master radial velocity template. Bottom: The cross-correlation peak 
obtained from the top spectrum using the master template. The Tonry-Davis R value for this 
star is 11.6, which is the median value of all the stars considered in the Phoenix radial velocity 
derivation 
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Fig. 6. — Radial velocities of Phoenix stars as a function of a number of parameters. Clockwise 
from top left: V r vs.V magnitude, V r vs.(V — I) color index, V r vs. position angle (in degrees), 
V r vs. galactocentric distance (in arcsec) within Phoenix. Crosses correspond to Phoenix set#l, 
open circles to Phoenix set #2.1 and solid circles to Phoenix set#2.2. Note that no trends in V r 
with respect any of these parameters, nor with the Phoenix set, exist. 
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Fig. 7. — Results of the ram-pressure stripping simulation: gray levels refer to the predicted HI 
velocity (darker corresponds to more negative values of the velocity) and white contours correspond 
to heliocentric radial velocities -27, -24, -21 and -18 km s _1 . The black contours show the 
predicted gas density (3.5, 3.0, 2.5, 2.0 and 1.5 xl0 19 cm~ 2 ) 



Fig. 8. — Top: Spectrum of a Carbon star in Phoenix (Phoenix set#2-6 in Table ||). Bottom: 
Spectrum of a young star in Phoenix (Phoenix set#2-10s3), probably a bright giant (luminosity 
class II) of early A spectral type. 
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Table 1: Observing log for the Phoenix and velocity standard stars 



Object a 2 ooo £2000 ARCHIVE FILE Exp Time (s) 



HD1 76047*1 


18' 


■59 


:44.12 


-34- 


28:14 


FOPS 


1999- 

. X U -J u 


-09- 


.14X23: 


:44: 


■30 


357 

•) 'J I 


fits 


10 


HD1 76047*2 

Lily ±_ 1 UUt: 1 -H—^i 


18: 


59 


:44.12 


-34: 


28:14 


FORS. 


1999- 


-09- 


14T23: 


:48: 


:29. 


001 


.fits 


20 


HD1 76D47 , », 
nu± i uui / sky#l 


1 8' 


59 


•45 7^ 


-34- 


98-14 


FOPS 


1 999- 


-D9- 


• 1 4T93 


CI 
. J- 


■ 0\J ■ 


499 


fits 

,11 Uo 


o\j\j 


HD176047#3 


18: 


:59: 


:44.12 


-34: 


:28:14 


FORS. 


1999- 


-09- 


15T00: 


:00: 


:21. 


.868. 


.fits 


150 


HD176047 sfcj/#2 


18: 


:59: 


:45.73 


-34: 


:28:14 


FORS. 


1999- 


-09- 


15T00: 


:06: 


:18. 


.654 


.fits 


600 


HD196983#1 


20: 


:41: 


:48.98 


-33: 


53:16 


FORS. 


1999- 


-09- 


15T00: 


:24 


:57. 


.582. 


.fits 


240 


HD196983 sfcj/#1 


20: 


:41: 


:47.44 


-33: 


53:16 


FORS. 


1999- 


-09- 


15T00: 


:34 


:39. 


793 


.fits 


300 


HD196983#2 


20: 


:41: 


:48.98 


-33: 


53:16 


FORS 


1999- 


-09- 


15T01: 


:21: 


:13. 


.492. 


.fits 


240 


HD196983#3 


20: 


:41: 


:48.98 


-33: 


53:16 


FORS. 


1999- 


-09- 


15T01: 


:28: 


:43. 


.540 


.fits 


120 


HD196983 sfcj/#2 


20: 


:41: 


:50.69 


-33: 


:52:57 


FORS 


1999- 


-09- 


■15T01: 


:34 


:41. 


.289 


.fits 


600 


HD203638 


21: 


:24: 


:07.99 


-20: 


:50:59 


FORS. 


1999- 


-09- 


15T03: 


:39 


:29. 


.409 


.fits 


5 


HD203638 sfcj/ 


21: 


:24: 


:09.27 


-20: 


:50:59 


FORS 


1999- 


-09- 


15T03: 


:42: 


:16. 


.512 


.fits 


600 


Phoenix set#l 


01: 


:51: 


:04 


-44: 


:26:05 


FORS. 


1999- 


-09- 


15T06: 


:03 


:47. 


.856 


.fits 


1500 


Phoenix set#l 


01: 


:51 


:04 


-44: 


26:05 


FORS 


1999- 


-09- 


15T06: 


:30 


:32 


.498 


.fits 


1500 


Phoenix set#l 


01: 


:51: 


:04 


-44: 


26:05 


FORS. 


1999- 


-09- 


15T06: 


:57: 


17. 


.324 


.fits 


1500 


Phoenix set #2.1 


01: 


:50: 


:57 


-44: 


:26:40 


FORS 


1999- 


-09- 


15T07: 


:34 


:53 


.935 


.fits 


1500 


Phoenix set#2.1 


01: 


:50: 


:57 


-44: 


26:40 


FORS. 


1999- 


-09- 


15T08: 


:01: 


:38. 


.649 


.fits 


900 


Phoenix set#2.2 


01: 


:50: 


:57 


-44: 


:26:40 


FORS. 


1999- 


-09- 


15T08: 


:18: 


:01. 


.062 


.fits 


1800 


Phoenix set#2.2 


01: 


:50: 


:57 


-44: 


:26:40 


FORS. 


1999- 


-09- 


15T08: 


:49: 


:45. 


.413 


.fits 


1800 


cpd432597#l 


06: 


:32 


:13.20 


-43: 


31:14 


FORS. 


1999- 


-09- 


15T09: 


:45 


:44. 


.645 


.fits 


12 


cpd432597#2 


06: 


:32: 


:13.20 


-43: 


31:14 


FORS. 


1999- 


-09- 


15T09: 


:48: 


:03. 


264 


.fits 


48 
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Table 2: Radial velocity standards. 



Star 


Vobs 


R 


a v 


Vhelio 


Vhelio - Vobs 


HD176047#1 


-40.8 


91.42 


4.9 


-42.8 


-2.0 


HD176047#2 


-46.1 


112.64 


4.0 


-42.8 


3.3 


HD176047#3 


-42.6 


115.85 


3.0 


-42.8 


0.2 


HD196983#1 


5.2 


62.45 


7.2 


-9.3 


-14.5 


HD196983#2 


-15.2 


98.44 


4.6 


-9.3 


5.9 


HD196983#3 


-10.1 


79.65 


5.6 


-9.3 


0.8 


HD203638 


10.6 


63.17 


7.1 


21.9 


11.3 



"The standard star radial velocities were taken from the following references: HD176047 and HD196983: Maurice et 
al. 1984; HD203638: Astronomical Almanac 1999 



Table 3: Dependence between wavelength zero-point shifts and telescope position. Values for the 
Phoenix sets are averages over all the slitlets 



Object 


Xobs 5577.34 - X obs 


Alt 


Azim 


HD176047 sfcj/#1 


5577.58 


-0.24 


79.6 


342.0 


HD176047 sfcy#2 


5577.45 


-0.11 


80.2 


358.8 


HD196983 sfcs/#1 


5577.70 


-0.36 


71.3 


304.0 


HD196983 sfcy#2 


5577.45 


-0.11 


80.2 


342.1 


HD203638 sfcy 


5577.15 


+0.19 


73.2 


99.4 


Phoenix set#l 


5577.66 


-0.32 


67.4 


334.4 


Phoenix set #2.1 


5577.05 


+0.29 


68.8 


18.8 


Phoenix set#2.2 


5577.08 


+0.26 


64.3 


34.6 
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Table 4: Summary of repeat velocity measurements 



Star 


Vhelio 


R 


< V > 


Vhelio- < v > 


HD176047#1 


-40.8 


91.4 


-43.2 


2.4 


HD176047#2 


-46.1 


112.6 


-43.2 


-2.9 


HD176047#3 


-42.6 


115.8 


-43.2 


0.5 


HD196983#1 


5.2 


62.4 


-6.7 


11.9 


HD196983#2 


-15.2 


98.4 


-6.7 


-8.5 


HD196983#3 


-10.2 


79.6 


-6.7 


-3.4 


Phoenix set#2.1-5 


-54.6 


8.8 


-66.3 


11.61 


Phoenix set#2.2-5 


-77.9 


13.4 


-66.3 


-11.61 


Phoenix set#2.1-8 


-26.3 


10.4 


-47.3 


21.01 


Phoenix set#2.2-8 


-68.3 


15.4 


-47.3 


-21.01 


Phoenix set#2.1-9 


-51.5 


8.1 


-52.3 


0.80 


Phoenix set#2.2-9 


-53.1 


21.1 


-52.3 


-0.80 


Phoenix set#2.1-12 
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Table 5: Data of the spectroscopically observed stars in Phoenix, with radial velocities for the RGB 
stars with sufficient S/N to obtain a well defined cross-correlation peak. 
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